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ABSTRACT 

The studies focused on the Thermally-Pulsing Asymptotic Giant Branch phase ex- 
perienced by low- and intermediate-mass stars are extremely important in many as- 
trophysical contexts. In particular, a detailed computation of their chemical yields 
is essential for several issues, ranging from the chemical evolution of galaxies, to the 
mechanisms behind the formation of globular clusters. Among all the uncertainties 
affecting the theoretical modelling of this phase, and described in the literature, it re- 
mains to be fully clarified which results are severely affected by the use of inadequate 
low-temperature opacities, that are in most cases calculated on the basis of the origi- 
nal chemical composition of the stars, and do not consider the changes in the surface 
chemistry due to the occurrence of the third dredge-up and hot-bottom burning. Our 
investigation is aimed at investigating this point. By means of full evolutionary mod- 
els including new set of molecular opacities computed specifically with the iESOPUS 
tool, we highlight which stellar models, among those present in the literature, need a 
substantial revision, mainly in relation to the predicted chemical yields. The interplay 
among convection, hot bottom burning and the low-temperature opacity treatment is 
also discussed. 
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1 INTRODUCTION 

After the end of helium burning in the core, stars with ini- 
tial masses M in the range 1 < M/Mq < 6 — 8 evolve 
through the Thermally-Pulsing Asymptotic Giant Branch 
(hereinafter TP-AGB) phase during which they experience 
a very rich nucleosynthesis whose products (e.g. He, C, N, 
Ne, Mg, s-process elements) may be convected to outermost 
layers and eventually ejected into the interstellar medium 
(ISM) by stellar winds, leaving a CO white dwarf remnant 
and an expanding planetary nebula. All these facts make 
the study of AGB stars extremely important for the chemi- 
cal evolution of the host system. 

In the last decades, the interest towards AGB evolution 
has increased essentially for two reasons: a) The observa- 
tions of high-redshift systems allowed the measurement of 
carbon and nitrogen abundances; b) AGB stars have been 
suggested as the main responsible for the self-enrichment of 
globular clusters, to explain the existence of multiple pop- 
ulations detected in almost all the clusters examined in the 
Galaxy on the basis of the photometric and s pectroscopic 
surveys ijNorris et alJll98ll : IVentura et al.ll200l[ ). 
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Chemical yields of AGB stars critically depend by the 
variations of their surface chemical composition, that may 
be produced not only by the first (and possibly the sec- 
ond) dredge-up, but also by two further mechanisms that 
are typical of the TP-AGB phase: a) The Hot Bottom Burn- 
ing (hereinafter HBB), i.e. the occurrence of proton capture 
reactions at the bottom of the external convective mantle, 
as a consequence of the high temperatures achieved by these 
layers in the most massive models; b) the Third Dredge-up 
(hereinafter TDU), by which the bottom of the envelope 
after each TP sinks inwards into layers previously touched 
by 3a nucleosynthesis. Both mechanisms alter the surface 
chemical composition: as for HBB we expect abundance 
changes resulting from the CNO-cycle reactions operating 
at equilibrium, whereas the TDU is associated to a net in- 
crease of the surface carbon and, to a smaller extent, of 
oxygen. 

The predictive power of the yields provided by the 
theoretical investigations is unfortunately undermined by 
the many uncertainties in the input-physics that affect 
the AGB description: mass loss, convection, treatment of 
the convective/radiative interfaces are unknown from first 
principles, and all have a strong impact on the physi- 
cal (and consequently chemical) evolution of AGB mod- 
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els jHerwidlioOoh . These problems are the reason for the 
great differences among t he results found by different re- 
search groups on this topic llMarigo 2001; Fc nner et al.l|2004l : 
IKarakas fc Lattanzioll2008l ; IVentura fc D'Antonall2005l V 

A further uncertainty, still to be fully investigated in 
this context, is the treatment of the molecular opacities, 
that are commonly calculated on the basis of the assumed 
metal content and the initial abundance distribution of the 
elements, neglecting any alteration of the surface chemistry 
associated to HBB and/or TDU. 

An e xploratory app roach related to this problem was 
made by iMarigcl l|2002h . who showed that when the C/O 
ratio exceeds unity, the formation of CN molecules deter- 
mines an increase of the opacity, that is expected to favour 
an expansion of the envelope, an increase of mass loss, and a 
faster consumption of the whole external mantle: this would 
diminish the number of TPs experienced by the star, and 
the effects of the TDU w ould be consid erably reduced. 

In a follow-up paper IMarigd (|2007T l considered the pos- 
sible effects driven by variations of molecular opacities in 
the most massive AGB stars (with M > 3.5M0), which 
experience both TDU and HBB. With the aid of envelope 
integrations it was found that if the dredge-up of carbon is 
efficient enough to lead to an early transition from C/O < 1 
to C/O > 1, then hot-bottom burning may be weakened, 
extinguished, or even prevented. 

A step forward in understanding the effects produced 
by the use of opacit y tables appropriate t o carbon-rich mix- 
tures was made bv lCristallo et ail (|2008l ). who presented a 
study focused on a 2Mq model of metallicity Z — 10 -4 . 
The authors concluded that the surface composition and 
the global yields by low-mass AGB stars of low metallicity 
are significantly affected by the use of the C-rich opacities. 
More extended sets of full TP-AGB evolutionary calcula- 
tions with variable molecular opacities have been recently 
presented by Weiss & Ferguson (2009) , covering wide ranges 
of stellar masses and metallicities. 

In the present paper we further extend the investigation 
on the effects of low-temperature opacities on the AGB evo- 
lution, so as to single out the regimes of stellar masses for 
which a major revision of the published models is required. 



2 DESCRIPTION OF THE MODELS 
2.1 The evolution code 

The models presented in this paper were calculated by 
means of the ATON code for st ellar evolution , a ful l de- 
scription of which can be found in I Ventura" et all l| 19981) . We 
provide here a brief summary of the numerical and phys- 
ical inputs used in the present computations. Convection 
was modelled according to the Full Spectrum of Turbu- 
lence (hereinafter FST) prescription bv lCanuto fc Mazzitellil 
|l99ll ). Mixing of chemicals and nuclear burning were cou- 
pled by mean s of a diffusive approach, following the scheme 
presented in ICloutman fc EolJ fl976); accordingly, extra- 
mixing was modelled by an exponential decay of convec- 
tive velocities beyond the formal borders, with an e-folding 
decay of I = C,Hp. During the two main phases of core nu- 
clear burning and in occurrence of the second dredge-up we 
used C, — 0.02, whereas in the TPs phase no extra-mixing 



was co nsidered. Mass loss was modelled according to lBlockerl 
1 19951 ). with the free parameter entering Reimers' prescrip- 
tion tjr = 0.0 2 . in a greement with the calibration given in 
IVentura et ail i^OPCD . 



2.2 Molecular opacities 

A key implementation in the present stellar models resides in 
the use of new tables of low-temperature opacities (1500 ^ 
T ^ 30 000 K), which are suitably constructed to follow the 
changes in the chemical composition of the envelope driven 
by TDU and HBB. 

A large set of tables of Rosseland mean (RM) opac- 
ities has been computed with a new tool, iESOPUS 
(Accurate Equation of State and OPacity Utility Soft- 
ware) described in Marigo & Aringer (2009; web-interface 
at http://stev.oapd.inaf.it/aesopus), to which the reader 
should refer to all the details. Each table of RM opac- 
ities covers a rectangular grid in the log(T)-log(i?) dia- 
gram, where the temperature is made vary in the interval 
3.2 log(T) 4.5, while the R = p(10 6 /T[K]) 3 parameter 
spans the range — 8 ^ log(ii) ^ 1. Following the formalism 
introduced by Marigo & Aringer (2009), the adopted ref- 
erence mixture assumes a metallicity Z Te { = 0.001, and an 
a-enhanced distribution of the elements, expressed by the 
parameter [a/Fe] = +0.4, with respect to the reference so- 
lar mi xture given by the compilation of IGrevesse fc Sauvail 
j 19981 ). 

According to the notation in Marigo & Aringer (2009), 
the assumed a-enhanced composition corresponds to mix- 
ture A, in which the abundances of a-elements (O, Ne, Mg, 
Si, S, Ca, and Ti) are enhanced while the metallicity is kept 
constant, Z — Z Ic f. The latter condition is fulfilled by requir- 
ing that the abundance variation of the enhanced elements 
is compensated by the total abundance variation of all the 
others. 

This implies that the concentrations of O and all other 
a-elements are incremented by ~ 23% compared to the 
scaled-solar values, while the abundances of the iron-group 
elements and all other metals, including C and N, are de- 
pressed by almost a factor of 2. As a consequence, the ref- 
erence chemical mixture - with Z ro f = 0.001, [a/Fe] = +0.4 
- corresponds to a carbon-to-oxygen ratio C/0~ 0.19, 
which is quite lower th an the (C/O)© ~ 0.49 of the 
IGrevesse fc Sauvail l|l998h solar mixture. 

As a next step, the reference a-enhanced mixture is 
further altered by varying the abundances of C, N and O, in 
order to account for the changes of these elements caused by 
TDU and HBB. To this aim we introduce three independent 
variation factors /, for C/O, C, and N defined as: 



v7 - fc/o 



Xc 
Xo 
X c 



Xc,, 



Xo, Te i 
fc Xc,rc{ 



where (Xc,ref, -^N,ref> Xo,ret) correspond to the abundances 
(in mass fraction) of carbon, nitrogen and oxygen in the 
reference chemical mixture with Z = Z la { = 0.001 and 
[a/Fe] = +0.4, while (Xc, Xs, Xo) denote the abundances 
of the same elements in the new mixture. By construction, 
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Figure 1. Rosseland mean opacity as a function of the temper- 
ature for several C/O values, as indicated. We assume log(ij) = 
—5. The abundance of carbon is made vary, while leaving un- 
changed that of oxygen. The reference chemical mixture is defined 
by Z Te{ = 0.001, X = 0.7, [a/Fe] = +0.4, and solar abundances 
from Grevesse & Sauval (1998). 

the variation factor of oxygen derives from condition: 
= foX ,rci with fo = fc/fc/o 

We notice that in this case the reference metallicity is not 
preserved, and the total metal content Z may be larger or 
lower than Z re f depending on the total CNO abundance vari- 
ation. 

The last composition parameter we need to specify is 
the hydrogen abundance X, and the normalization condi- 
tion is fulfilled by requiring that the helium abundance is 
given by Y = 1 — X — Z. The assumed values of the four 
independent composition parameters X, fc/o, fc and fo 
are listed in Table [1] In total, 1296 opacity tables have been 
computed for all combinations of these parameters. For its 
relevance to the opacity results, the last column of Table [1] 
shows also the C/O ratio corresponding to each fc/o- 

In Fig. [T] each curve shows the predicted RM opaci- 
ties for a given value of the C/O ratio, spanning over three 
orders of magnitude, passing from ~ 0.006 to ~ 6. While 
for log(T) > 3.55 the different curves almost overlap as 
the Rosseland mean is controlled by the continuum opac- 
ity of hydrogen (bound- free and free- free transitions), for 
log(T) < 3.5 they start to deviate significantly. As discussed 
in Marigo & Aringer (2009), for 3.55 < log(T) < 3.45 the 
RM opacity is mainly affected by the anion H~, Thom- 
son electron scattering, and CN molecular bands, while for 
3.45 < log(T) < 3.2 the most important absorbers are 
the oxygen-bearing molecules such as H2O, VO, TiO for 
C/O < 1, while C-bearing molecules such as CN, C2, C3, 
HCN, and C 2 H 2 control the opacity for C/O > 1. 

The results presented in Fig. [T] provide a case overview 
of the changes in opacities driven by changes in C/O con- 



Table 1. Values of the composition parameters fi adopted in 
the computation of the RM opacity tables. The sequence of the 
variation factor fc/o i s designed to allow a good sampling of the 
critical points of C/O at which significant opacity changes are 
expected. Refer to the text for more explanation. 



X log(/ c/ o) log(/c) M/n) C/O 



0.5 


-1.50 


-1.5 


0.0 


6.16610 


0.7 


-1.00 


-0.5 


0.5 


1.950 10 


0.8 


-0.50 


0.0 


1.0 


6.16610 




0.00 


0.5 


1.5 


1.950 10 




0.30 


1.0 


2.0 


3.890 10 




0.50 


1.5 


2.5 


6.16610 




0.70 






9.77210 




0.75 






1.096 




0.80 






1.230 




1.00 






1.950 




1.30 






3.890 




1.50 






6.166 



sequent to TDU and HBB. Assuming C/O rc f = 0.19 as the 
reference value, the opacity curves with C/O < 0.19 cor- 
respond to an effective decrease of the carbon abundance 
(log(/ c ) < 0), hence representing the cases in which HBB 
burns carbon in favor of nitrogen; on the other hand the 
opacity curves with C/O > 0.19 correspond to an effec- 
tive increase of the carbon abundance (log(/ c ) > 0), hence 
suitable to describe the cases in which TDU prevails. 

One striking feature concerns the different sensitiveness 
of the RM opacities to changes in carbon abundance. The 
diversification among the opacity curves is notably larger 
in mixtures with C/0> 1, while it keeps rather small as 
long as C/0< 0.9, at least for log(T) ~ 3.4 - 3.6. As dis- 
cussed in Sect. [3] this feature has important consequences 
for the behaviour of the TP-AGB models, i.e. C-rich models 
are in general more importantly affected by changes in the 
chemistry than O-rich models. Finally, we remark that the 
adopted set of fc/o and fc accounts for both a decrease and 
an increase of C/O and C (hence also O). Our choice has a 
well-founded motivation. In fact, during the TP-AGB evolu- 
tion newly-synthesized C and O may be not only convected 
to the surface by TDU, but their envelope abundances may 
be also depleted by HBB, i.e. partially burnt in favour of N 
by the CNO-cycle. In this case the evolution of the surface 
C/O ratio does not simply trace the increase of C due to 
TDU, being rather the result of concurring and competing 
effects (TDU and HBB), which should be considered in RM 
opacities suitable for AGB modelling. 



3 THE ROLE OF THE NEW OPACITIES 

The evolutionary sequences presented in this work have 
been followed from the pre-MS through the whole TP-AGB 
phase, up to the latest stages of the AGB evolution, when 
almost all the envelope was lost. The initial chemical com- 
position of the models is typical of an intermediate metal- 
licity globular cluster: Z=0.001 and Y=0.24, with an alpha- 
enhanced Qa/Fe]=0.4) distribution of the elements; the in- 
dividual abundances of the various species are taken from 
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Table 2. Chemical yields and evolutionary properties of AGB stars 



jESOPUS low-T opacities 
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Figure 2. The evolution of the surface C/O ratio in three stellar models calculated with the standard opacities (solid line), and with 
the new opacities that account for the changes in the CNO abundances (dotted line). 



iGrevesse fc Sauval (|l998l ). We focus our attention on three 
models of mass M = 2.5, 3.0, 3.5 Mq, as their evolution en- 
compass the whole range of cases that are relevant for the 
present investigation. Each evolutionary model was calcu- 
lated with the iESOPUS opacities assuming both the refer- 
ence chemical mixture, hence neglecting the surface chem- 
istry changes during the TP- AGB phase ("standard case"), 
and accounting for the CNO abundance variations ("new 
case"). 

The chemical content of the ejecta, averaged over the 
whole AGB phase, are indicated in Table 2; R(CNO) in col.6 
is the ratio between the overall CNO content of the ejecta 
compared to the initial value. Table 2 also includes, in the 
last two columns, the maximum temperature reached by the 
bottom of the convective envelope ILce an d the radius of the 
star before the envelope is lost. We limit our discussion to the 
CNO elements here, because this allows to better disentangle 
the effects of the opacity description from those associated to 
the cro ss sections of the p roton-capture reactions by heavier 
nuclei l|lzzard et al.ll2007l ). 

Figure [2] shows the evolution of the surface C/O ratio 
of the models. The 2.5M0 model evolves as a carbon star. 
The use of the iESOPUS CNO opacities favours larger radii 
(see last column of Table 1), with a consequent increase 
of the mass loss rate: the faster consumption of the enve- 
lope prevents the C/O from exceeding 2, at odds with the 
model calculated with the standard opacities, whose C/O 
ratio approaches 3 in the final stages of the evolution. The 
average C/O of the ejecta is consequently ~ 40% higher in 
the standard case. The cooling effect of the CNO opacities is 
confirmed by the differences in the maximum temperature 



attained by the bottom of the envelope; note that these tem- 
peratures are well below the threshold (~ 60 — 70K) neces- 
sary to activate HBB, and this is the reason why the nitrogen 
content of the ejecta is the same in the two cases. 

The 3.OM0 shows up the most striking differences (see 
the middle panel of fig[2} . In the standard case we see an ini- 
tial increase of the surface C/O, followed by a phase of car- 
bon depletion, that is the signature of HBB. The evolution 
of C/O in the jESOPUS CNO model is completely different: 
the increase of the surface carbon stops when C/0~ 1.5, 
and it keeps approximately constant for the remaining evo- 
lution, the small depletion during the CNO burning phase 
being compensated by the increase determined by the TDU 
episodes following each TP. In this case the cooling of the 
envelope (see the different temperatures of the two models 
in the last column of Table 1) prevents HBB, so that this 
latter is extinguished. Unlike the lowest masses case, the 
yields of the standard model have a lower C/O ratio; also 
the nitrogen content of the ejecta is different, because HBB 
favours a larger increase of the nitrogen abundance in the 
standard case. 

The right panel of fig[5]shows that in the M = 3.5M@ model 
the differences introduced by the CNO opacities are small. 
This is motivated by the low C/O ratio, never exceeding 0.5. 
The decrease of C/O after a few TPs indicates that HBB is 
active in both models. These results suggest that HBB, when 
strong, is efficient in both cases. This fact is explained con- 
sidering that in the range 0.05 < C/O < 0.5, spanned by the 
M = 3.5M TP-AGB model, the predicted RM opacities 
differ relatively little, as displayed by the curves in Fig. [1] 
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4 THE DEMAND FOR A REVISION OF THE 
AGB MODELLING 

The results presented in the previous section show that us- 
ing the appropriate low temperature opacities, that accounts 
for the formation of new molecules when the C/O ratio ap- 
proaches and exceeds unity, may change substantially the 
evolution of the surface abundances of carbon, nitrogen and 
oxygen during the AGB evolution. In low-mass AGB stars, 
when HBB is absent, the use of fixed opacities largely over- 
estimate the carbon yield and the C/O ratio. On the other 
hand, we find that for the models achieving mild HBB con- 
ditions, the treatment of the opacities is relevant to main- 
tain or quench HBB itself, so that, unlike the previous case, 
we expect much higher yields of 12 C, and smaller 14 N abun- 
dances. This effects is restricted to a limited range of masses, 
~ 0.5Afo, because the results for more massive models in- 
dicate that when a stronger HBB is present a depletion of 
the surface carbon is achieved independently of the opacity 
treatment; these models show some differences in the oxygen 
yields (that is depleted more efficiently in the standard case) , 
but even this latter difference is expected to vanish with in- 
creasing mass. These results are in full agreement with the 
findings by Marigo (2007), obtained by means of synthetic 
TP-AGB calculations based on numerical integrations of 
complete envelope models, which included an approximate 
description of variable molecular opacities (Marigo 2002). 

The question of which TP-AGB models, among the 
most massive ones, require substantial revision, as a conse- 
quence of introducing adequate molecular opacities, is there- 
fore shifted to the well-known and long-standing problems 
associated to the AGB evolution: which conditions are nec- 
essary to achieve HBB l|Ventura fc D'Antonall2005l '). It has 
been shown that HBB is strongly connected with the treat- 
ment of convection, and that when this latter is modelled 
efficiently (e.g. when the FST model i s adopted) HBB is 
favoured l|Blocker fe SchoenberneJll99ll ). The results of the 
previous section indicate that the range of masses whose 
yields are crucially affected by the new opacities differs de- 
pending on the convective modelling: when this latter is effi- 
cient, only models for masses smaller that ~ 3 — 3.5M© the 
use of the appropriate opacities is mandatory not to substan- 
tially mistake the estimation of the yields, whereas AGB 
models computed with a low-efficiency convective model 
would demand even a larger revision, extended to almost 
all the intermediate masses. 



5 CONCLUSIONS 

In this paper we focused on the role played by the treatment 
of low-temperature opacities, consistently coupled with the 
surface chemical composition, on the description of the AGB 
evolution. The scope was to clarify what part of the existing 
literature, and commonly used in many astrophysical con- 
texts, needs a severe revision due to the frequent incorrect 
estimate of the surface opacities. To this scope, we present 
new evolutionary sequences at low metallicity (Z = 0.001), 
that cover most of the situations typical of the AGB phase 
in terms of TDU and HBB, and compare the standard re- 
sults with the new models, calculated with a set of opacities 
appositely built for the present investigation, that accounts 



for the variation of the surface chemistry. Our main findings 
are the following: 

(i) In absence of HBB, the use of correct opacities is 
mandatory, otherwise the carbon yields and the C/O ra- 
tio of the ejecta are largely overestimated. This holds for all 
the low masses, that achieve TDU but not HBB. 

(ii) For a narrow range of masses, those achieving only 
mild HBB conditions, the necessity of using the new opac- 
ities is even more urgent, because the surface evolution de- 
pends dramatically on the opacity treatment: in this case 
HBB is quenched, carbon is never depleted if not in mi- 
nor quantities during the interpulse phase, so that the tra- 
ditional computations strongly underestimates the carbon 
yields. 

(iii) For more massive models, strong HBB prevents the 
surface C/O ratio to increase substantially, so that the opac- 
ity treatment seems not critical for the results obtained. 

The general validity of these conclusions should be further 
tested by extending the calculations at both higher and 
lower metallicities. This work is in progress. In any case, 
thanks to the present availability of RM opacity tables for 
arbitrary chemical mixtures, (i.e. Marigo & Aringer 2009 
(the iESOPUS web-tool), Helling & Lucas 2009, Lederer & 
Aringer 2009), the use of molecular opacities, consistently 
connected to the actual surface chemistry, could become the 
standard choice for stellar evolutionary calculations to come. 
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